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Abstract. In calorimetric neutrino mass experiments, where the shape of a beta decay spectrum has to 
be precisely measured, the understanding of the detector response function is a fundamental issue. In 
the MIBETA neutrino mass experiment, the X-ray lines measured with external sources did not have 
Gaussian shapes, but exhibited a pronounced shoulder towards lower energies. If this shoulder were a 
general feature of the detector response function, it would distort the beta decay spectrum and thus mimic 
a non-zero neutrino mass. An investigation was performed to understand the origin of the shoulder and its 
potential influence on the beta spectrum. First, the peaks were fitted with an analytic function in order 
to determine quantitatively the amount of events contributing to the shoulder, also depending on the 
energy of the calibration X-rays. In a second step, Montecarlo simulations were performed to reproduce 
the experimental spectrum and to understand the origin of its shape. We conclude that at least part of 
the observed shoulder can be attributed to a surface effect. 
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1 Introduction 

Low temperature microcalorimeters are expected to play 
a key role in future direct neutrino mass measurements 
PQ. With these detectors it is possible to perform very 
sensitive beta end-point studies in a calorimetric configu- 
ration, i.e. with the beta source embedded in the detector. 
In this way, the detector measures all the energy released 
in the beta decay except that carried away by the neutrino, 
and many systematic effects showing up in other configu- 
rations are avoided. Nevertheless past neutrino mass ex- 
periments with low temperature microcalorimeters have 
shown how critical is the understanding of the detector 
response function. 

In past years few calorimetric experiments using 187 Re 
as beta decaying isotope have been carried out both in 
Milano (MIBETA) and in Genova (MANU). A final sensi- 
tivity on the neutrino mass of around 15eV was achieved 
by both groups [2U3B1E]- in these experiments the neu- 
trino mass squared to 2 is measured by looking at the tiny 
deficit of events in an energy interval below the beta decay 
end-point as small as few times the neutrino mass itself. 
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For this purpose, the experimentally measured beta decay 
spectrum has to be interpolated by the function obtained 
as convolution of the theoretical beta decay spectrum with 
the detector response function F(E,E ) - defined as the 
measured response for the deposition of a fixed amount 
of energy Eq. A precise evaluation of F(E,E ) is there- 
fore mandatory to avoid systematic uncertainties. The re- 
sponse function F(E,Eq) is expected to be a symmetric 
Gaussian line centered at Eq whose width is caused by 
the electronic noise and is independent of Eq. In practice, 
however, the width of the response function depends on Eq 
and is larger than expected assuming only the electronic 
noise contribution. Moreover, often the line is not sym- 
metric and, even worse, its shape may depend on where 
Eq is deposited inside the detector. Usually the response 
function F(E,Eq) is investigated by means of suitable 
calibration sources, but this may be a source of system- 
atic effects. First of all, it is almost impossible to have 
really mono-energetic sources: the understanding of the 
calibration source emission spectrum is therefore part of 
the game. Second, it is very difficult to emulate with a 
calibration source the beta decay energy deposition: the 
calibration source has to cover the same energy range as 
the beta particles, the energy depositing interactions must 
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be spatially equivalent both in their global distribution 
throughout the detector and in the single tracks, and the 
energy deposition must be done by the same type of parti- 
cle, i.e. an electron. Last but not least, it must be possible 
to remove the calibration source during the beta decay 
measurement without affecting the detector. 

In calorimetric neutrino mass experiments, given the 
low energy end-point of the beta decays (e.g. 187 Re end- 
point is 2.465 keV 3 ), the only viable approach is to 
use an external low-energy X-ray sourc^E The energy de- 
positions are then caused mostly by the primary photo- 
electrons and by the cascade of secondary X-rays and 
Auger electrons. Few differences to the beta decay electron 
interactions are immediately apparent. First, whereas the 
creation of a photo-electron by an X-ray is followed by 
a cascade of secondary atomic radiation (mostly Auger 
electrons) , low energy betas produce very little secondary 
radiation (mostly photons). Second, with the X-rays it is 
experimentally difficult to achieve a uniform illumination 
of the active detector volume as it happens with the beta 
particles. It is therefore extremely important to under- 
stand the measured response function F(E,E ) in order 
to disentangle the contributions to its shape caused by the 
X-rays from an external source. 

The first aim of the present work is to obtain a precise 
functional description of the detector response function 
F(E,Eq) measured for X-ray absorption. Then we try to 
establish whether and how this response function can be 
used to analyze the beta decay spectrum. To accomplish 
this, the physical processes underlying the response func- 
tion for X-ray absorption will be investigated comparing 
different hypotheses with the experimental data. Further- 
more, processes taking place in the source assembly that 
can cause low energy tails have to be carefully investi- 
gated. Once a model is confirmed, Montecarlo simulations 
can be used to indirectly obtain the energy response for 
the beta decay electrons. 



2 Experimental Set-up and Data 

The present work is based on the data collected during the 
activity of the MIBETA experiment from 2000 to 2004 [1 
[3]. In particular, the three runs named RUN9, RUN14 and 
RUN15 are considered. These three runs were carried out 
with the same detectors but with different configurations 
because of their different goals (see Table Q] for details 
on the set-ups). The detectors are ten microcalorimeters 
made of small AgReC>4 crystals glued to doped silicon chip 
thermometers. The crystal mass ranges from 250 to 300 
fig, for a total mass of about 2.68 mg (see [3] for details on 
the set-up). The major differences between the three runs 
are their durations and the calibration source assemblies, 
although they share the same low energy X-ray calibra- 
tion source based on the fluorescence of low Z materials 
exposed to the X-rays of two 55 Fe sources. 

1 The use of an external low energy monochromatic beta 
source would be much more problematic as shown by spectro- 
metric neutrino mass experiments. 



RUN Calibration Calibration sources Usable 

ID time [h] 55 Fe Pb shield 44 Ti detectors 

RUN9 240 yes no no 8 

RUN 14 1000 yes yes no 8 

RUN 15 1000 yes yes yes 2 

Table 1. Summary of the experimental conditions in the three 
runs considered in this work. The 55 Fe calibration source refers 
to the fluorescence source described in the text. The last col- 
umn reports the number of detectors which were used for the 
analysis presented in this work (see text). 



RUN9 was intended to be the MIBETA high statistics 
measurement with the goal of achieving the best statisti- 
cal sensitivity on the neutrino mass and it features a first 
unshielded version of the fluorescence calibration source. 
The measurement was stopped after only 1000 hours be- 
cause of the high continuum background caused by the 
55 Fe inner Bremsstrahlung. In RUN14 a lead shielding 
for the 55 Fe sources brought a drastic reduction of the 
background and the high statistics measurement was com- 
pleted as planned. RUN15 was carried out with an addi- 
tional 44 Ti source with the aim of investigating the detec- 
tor energy response. 

The fluorescence source is made of two primary 5 mCi 
55 Fe sources irradiating two composite targets containing 
Al, CaF2, Ti, and NaCl. Therefore, the detectors are ex- 
posed to the K fluorescence lines of Al, Ca, Ti and CI as 
well as to the Rayleigh scattered K X-rays of Mn result- 
ing from the electron capture decay of 55 Fe. Thanks to a 
mechanism operated from outside the cryostat, the 55 Fe 
sources can be shut off in order to stop the emission of the 
fluorescent X-rays and to measure the beta decay spec- 
trum without background. In RUN9 the 55 Fe sources are 
shut off by two copper foils, while in RUN14 and RUN15 
the 55 Fe sources are moved inside a thick lead shielding. 
In all runs, additional weak fluorescence lines are observed 
because of the several materials present in the set-up, like 
the Pb used for the primary source housing in RUN14 and 
RUN15 or the stainless steel used for the primary source 
rails and case. In particular, in the RUN14 spectrum - 
shown in Figure [T]- the Pb M lines below 2.5 keV and the 
Cr Ka line at 5.415 keV are evident. Table [2] lists the rel- 
evant X-ray energies together with their respective ranges 
in AgReC-4. 

The data of RUN14 used for the present work are the 
subset collected during the calibration cycle of the mea- 
surements with only eight of the ten detectors. The total 
measurement time amounts to about 1000 hours and the 
resulting spectrum is shown in Figure [T] The beta spec- 
trum of 187 Re has been subtracted - about 6 x 10 5 counts 
between 1.0 and 2.465 keV- to enhance the Al lines and 
ease their investigation. Unfortunately their shape has 
been altered by the software cuts applied to the data in 
order to separate spurious events from true beta decays 
(see [5] for details about the data analysis). It can be 
clearly seen in Figure Q] that all lines deviate from Gaus- 
sian shapes, clearly exhibiting a tail towards lower ener- 
gies. In spite of the lower statistics, Figured] shows that 
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Fig. 1. Summed calibration spectrum of eight detectors from 
the RUN14 measurement. The beta spectrum has been sub- 
tracted (see text). 




energy [keV] 



Fig. 2. Ti Ka lines from the calibration spectrum of each 
channel which is included in the sum of Figure[T] In spite of 
the lower statistics all spectra show the low energy tail observed 
in the summed spectrum. 



also the lines in the individual channel spectra have the 
same low energy tail. The FWHM energy resolution on the 
CI K Q peak - the closest to the beta decay end-point - is 
about 29 eV. The same tails are actually observed every 
time calibration spectra with high statistics are collected 
with AgRe04 detectors. 

As a first step to understand the observed calibration 
spectrum, two Montecarlo simulations were performed us- 
ing the Geant4 toolkit [7J. The first simulates simply the 
direct interaction of mono-energetic X-rays in the AgReC>4 
absorbers showing that no tails are expected in the fluo- 




Fig. 3. Virtualized representation of the experimental setup 
with the Fe source and the fluorescence target. 
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Ti K/3 
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2.6 


Ti Ka 
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4.014 
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Ca Ka 
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CI K/3 


2.818 


1.0 


CI Ka 
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Pb M/3 
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Table 2. Energies ,6] and ranges of calibration 7 and X-rays 
(only Kai and K/3 are listed). 





peak position 


Rhenium X- 
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line name 




energy 




range 




[keV] 


[keV] 


[Mm- 1 
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Re Ka 1 escape 


17.196 


61.141 


0.044 


23 


Re Ka 2 escape 


18.618 


59.717 


0.035 


28 


Re K/3 1 escape 


9.027 


69.31 


0.024 


42 



Table 3. Escape peaks: positions together with the energies 
[5] and ranges of X-rays involved in the escape process. 




energy [keV] 

Fig. 4. This simulation of the Fe fluorescence spectrum 
shows that no shoulders are expected for the fluorescence 
peaks. The shoulders in the Mn peaks stem from Compton 
back-scattering. 

rescence peaks due to mechanisms like emission and loss 
of secondary radiation following the photoelectric absorp- 
tion. 

The second is a detailed Montecarlo simulation of the 
whole experimental set-up. Figure[3] shows a virtualized 
representation of the simulated set-up. The spectrum in 
Figure^] was obtained using the Penelope low energy ex- 
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tension of Geant4 [7]. The standard Geant4 low energy 
extension was discarded since it introduced in the spec- 
trum extra features which were not seen in the experimen- 
tal dat£0- The Montecarlo spectrum corresponds to about 
2 x 10 11 generated 55 Fe decays - i.e. to about one hour 
of experimental live time - including all X-rays and all 
Auger electrons. The spectrum in Figure|4] clearly shows 
that no tails are expected in the fluorescence peaks due 
to scattering in the set-up. At the same time the simu- 
lation shows that a large exponential tail is expected for 
the Mn Rayleigh scattered peaks due to a small fraction 
of Compton back-scattering events. The exponential tail 
contains about 10% of the total Mn Ka events and has a 
decay constant A of about 10keV _1 . 

In order to understand the origin of the tail observed 
in all other peaks, as mentioned above, it is important 
to take into account the different spatial distributions of 
X-ray and 187 Re decay events. As can be seen from Ta- 
ble[H the calibration X-rays have ranges in AgRe04 of 1-5 
/j,m only. As the AgReC>4 crystals have an average thick- 
ness of 240 ^m, this means that the energy deposition of 
the X-rays happens only in a thin surface layer containing 
less than 2% of the total crystal volume, while the beta 
decays occur uniformly across the whole crystal volume. 
The purpose of RUN15 was to experimentally investigate 
whether such a tail occurs also for events depositing en- 
ergy inside the crystal volume like the 187 Re beta decays. 
To accomplish this, a 44 Ti source was introduced into the 
set-up. 44 Ti decays by electron capture to 44 Sc emitting 
also one 7-ray with an energy of about 78.337 keV and 
a branching ratio of about 98%. Its energy is therefore 
just above the K-edge of Rhenium at 71.7keV, so peaks 
at about 18, 17 and 9keV due to the escape of Kai, Kc*2 
and K/3 Re X-rays, respectively, are expected in the exper- 
imental spectrum. Since for an energy of about 78 keV the 
attenuation length in AgReC>4 is larger than 200 /im the 
escape peaks, in spite of being located in the low energy 
region of the spectrum, are the results of interactions dis- 
tributed almost uniformly in the crystal. The energies and 
ranges related to the escape process are listed in Tabled 
To excite the escape peaks a 7-ray is much preferred over 
a X-ray in order to limit the peak broadening to just the 
contribution of the Re atomic transition natural width. 
The spectrum in Figure[5] is a Montecarlo simulation of 
the exposure to a 44 Ti source of a AgReC>4 detector with 
size similar to the ones considered here. It is possible to 
recognize the peaks due to the escape of Kai, Kev2 and K/3 
Re X-rays mentioned above: the escape peaks between 5 
and lOkeV reflect the complex structure of the K/3 X-ray 
line. The flat continuum is caused by Compton scattering 
of the 44 Ti 7-rays. 

In RUN15 the 44 Ti source could not be shielded and 
the detectors were always exposed to its radiation. This 
caused a high level of background during the measure- 
ments with the fluorescence source (see Figure More- 

2 In order to obtain correct results from Geant4, we have 
modified the Livermore Evaluated Atom Data Library (EADL) 
fluorescence files for the atomic number of interest because the 
original one had wrong X-ray transition energies. 



10 % N 



Ka escape : 




energy [keV] 

Fig. 5. Simulation of the 44 Ti spectrum showing the escape 
lines. 




energy [keV] 

Fig. 6. Summed spectrum of channels 5 and 6 from the 
RUN15 44 Ti measurement with the Re escape peaks. The beta 
spectrum has been subtracted - about 1.5 x 10 counts between 
1.0 and 2.465 keV. 



over, because of the intensity of the Ti source the anal- 
ysis of RUN15 data is severely impaired by pile-up. Due 
to the asymmetrical placement of the 44 Ti source, two of 
the eight detectors (detector 5 and 6) show a lower pile-up 
rate with a smaller impact on their performance. There- 
fore, for the analysis discussed in this paper only these two 
detectors were used. Figure H2 shows the energy spectrum 
obtained from the whole RUN15 measurement for these 
two detectors. The total measurement time is again about 
1000 hours. The FWHM energy resolution is about 31 eV 
for the CI K a line, which is only 7% worse than in RUN14. 

As discussed in the following, to investigate the Al and 
CI peaks we have used also data from RUN9. The total cal- 
ibration time in this case amounts only to about 220 hours. 
Along with a lower statistics, the RUN9 spectrum has a 
slightly worse energy resolution with respect to RUN14. 
On the other hand, RUN9 is free from the low energy X- 
ray peaks from Pb fluorescence and gentler offline cuts did 
not affect the shape of the Al peak. 

In the following we discuss the fitting of the experimen- 
tal calibration spectra as well as the spectrum of the 44 Ti 
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measurements. Then we present the results of the peak 
shape investigation by means of Montecarlo simulations. 



3 Fitting of the Experimental Spectra 

Tails in experimental spectra have also been observed in 
standard ionization semiconductor detectors. In particu- 
lar, in alpha silicon detectors a tail is attributed to energy- 
loss fluctuations in the collisions of the alpha particles 
with the silicon nuclei in the crystal [8]. The resulting 
asymmetric peaks T(E, E ) are generally described as the 
convolution of a normalized Gaussian G(E,Eq) with an 
exponential R(E, E ) 



G(E,E ) = 



1 



exp 



1 ( E-Ep 

2 V a 



CTV27T 

R{E, E Q ) = A exp [{E - E )X] ■ h{E - E) 
T{E,E ) = G(E,E )®R(E,E ) 

2' 



(1) 

(2) 
(3) 



A 

2 exp 



(r-r„,A- 



•erfc 



E — En 



aX 



aV2 y/2. 



where Eq is the peak position, a is the Gaussian standard 
deviation, A is the exponential tailing constant, h(Eo — E) 
is the Heavyside function, and erfc is the complementary 
error function. 

The above function T(E, E ) describes one Gaussian 
peak with an exponential tail. However, we found that 
this simple approach does not provide a satisfactory de- 
scription of the lines in our spectra (Figure [T]). In our 
case, the best description is given by the sum of three 
independent components: one main symmetric Gaussian 
peak G(E,Eq) and two asymmetric Gaussian peaks like 
the above T(E,Eq). The two asymmetric peaks are the 
convolution of the Gaussian G(E,Eq) with two exponen- 
tials, Ri(E,Eq) and R 2 {E,Eq), with tail parameters Ai 
and A2: 

F(E, E ) = G(E, E ) <8> [6(E) + R 1 (E, E ) + R 2 {E, E )} 

2" 



F(E,E ) = ^exp 



1 ( E-E p 



(4) 



Aexpi y exp 



•erfc 



(E — E )Xi 
E — Eo <j Ai 



A 2 

+ A exp2 — exp 



ay/2 y/2_ 

(E — E )X2 



aXi 

7f 



cA 2 
72 



•erfc 



E — Eo & A2 



aV2 y/2 



Acauss is the area of the Gaussian peak, A exp i and A exp 2 
describe the amplitudes of the two exponential tails. As 



will be argued in the following, one of the two exponential 
tails can be explained as a consequence of the shallow X- 
ray absorption depth. 

In spite of the relatively large width of the Gaussian 
peaks (ranging from about 25 to about 45 eV FWHM) a 
precise description of the measured spectra calls for a de- 
tailed modeling of the X-ray emission. This is particularly 
true for the Re Ka escape peaks whose natural line width 
r - about 47 eV - is comparable to the Gaussian line 
width a - about 28 eV. Therefore, both the doublet struc- 
ture of the Ka lines and the natural widths r of X-ray 
transitions have to be included. 

To take into account the natural line widths, for each 
peak the above response function F(E,Eo) was numer- 
ically convoluted with an asymmetric Breit-Wigner for- 
mula [TU] : 



BW(E,E ,r,5 AS ) = 
for E < E 

BW(E,E ,r,S AS ) = 
for E > E 

r R = 



Sas 



El 



tt(1+^s) (E-E ) 2 +El/A 



(5) 



En 



tt(1+5as) (E-Eo) 2 +E 2 R /A 



2r 



1 + S 



r L = 2r- r E 



AS 



(6) 
(7) 



Hereby, r is the natural line width and 8 as = -Tl / Er the 
asymmetry index. The values used for the fits are listed in 
TableUJ It is worth noting that it is not quite straightfor- 
ward to find reliable atomic parameters in the literature, 
especially for the low Z atoms like Al and CI. Moreover, 
line widths and asymmetry indexes depend on the chem- 
ical form of materials [PTIH^] . Whenever possible we have 
used the most recent and appropriate experimentally de- 
termined parameters. Lacking a direct experimental de- 
termination, the natural widths of CI and Re lines are the 
ones recommended in a recent compilation of experimental 
and theoretical values. When published values were miss- 
ing - as for Al, CI and Re - the asymmetry index has been 
arbitrarily set to 1.0. Given all the above, some systematic 
uncertainty on the fit parameters has to be expected. 

The Re Kal and Ka2 escape lines and the K/3 lines - 
whose structure is too fine and complex to be appreciable - 
can be satisfactorily fitted with the single peak H S (E, Eo) 
fit function: 

H S (E,E ) = F(E,Eo)®BW(E,E ,r,S AS )+background 

(8) 

For the other Ka lines, the Kal-Ka2 structure - although 
not fully resolved by our detectors - must be accounted 
for by the doublet fit function Hd(E,Eo): 

H d {E,E ) = F{E,E )®[BW(E,E Kal ,r Kal ,6 A s,K a i) 

BW{E,E Kal - AE,r Ka2 ,S AS , Ka2 )} 



LKal 



background 



(9) 



The energy difference AE and relative intensities of the 
Ka lines Ixai and Iku2 are listed in Tabled as well. Since 



6 



E.Ferri et al.: Eur. Phys. J. A (2012) 48:131 



tails from peaks at higher energies add up to an approxi- 
mately constant background at lower energies, a constant 
background parameter was added in both H s (E,Eq) and 
Hd{E,E ). The single and doublet peaks of the summed 
spectra in Figure[TJwere fitted individually using H s (E, E ) 
and Hd(E,E ), respectively. The free fit parameters are 
the peak position Eq (or Eku\ for doublets), the three 
amplitudes AQ auss , A exp \ and A exp 2, the Gaussian stan- 
dard deviation a, the exponential tail parameters Ai and 
A2 and the constant background level. Because of the low 
statistics and the large number of fit parameters, the K/3 
peak fit results have not been considered. 

Figure [7] shows as examples the fits of the Ti Ka and 
the Mn Ka peaks for the RUN14 summed spectrum. The 
bump on the right side of the Mn peak is caused by the 
K/3 peak of Cr at 5.950 keV, also included in the fit. 

FigurejEHH] andHU] show the results for the parameters 
a, Ai, A2, A cxp i and A cxp 2 as found by fitting the Ka 
peaks in the calibration spectra. In these plots the tail 
amplitudes A exp \ and A exp 2 are normalized to the total 
peak area Ai Q t — AQ auss -t- A exp \ -t- A exp 2. 

As far as the tail parameters are concerned, Figure[5] 
and [10] show the fit results only for a subset of the peaks 
we have analyzed. These results were selected considering 
the stability of the fit procedure with respect to the initial 
parameter values and to the energy interval as well as the 
convergence to meaningful values of all free parameters. 

According to this considerations, a clean analysis of 
RUN14 data is possible only for the Ca, Ti and Mn peaks. 
As for the CI peak, in the RUN14 spectrum the two M 
Pb peaks are too close for a correct estimation of the 
tail parameters (see spectrum in Figurc[JJ), so the RUN9 
spectrum was used instead. However, the low statistics 
and high background in this measurement prevent from 
properly estimating the long tail parameters, which are 
accordingly omitted in Figure|H] and [TO] For the Al peak 
in RUN14, the determination of the long tail parameters 
is not possible because the peak shape has been altered 
by the software cuts applied to the data in order to sep- 
arate spurious events from true beta decays (see [3] for 
details about the data analysis) . Although in the analysis 
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Fig. 11. Comparison of the normalized Ti Ka peaks in RUN14 
(black) and RUN15 (grey). 



of RUN9 gentler offline cuts did not affect the shape of the 
Al peak, the long tail parameters are negatively influenced 
by the presence of the beta spectrum. Therefore, they are 
also omitted in Figurel9l and ITOl 

The results of the fits for the a parameter are dis- 
played in Figure[SJ The a parameter is the only one show- 
ing an evident trend with energy. Rs energy dependency 
may be described by a sum of three squared contributions 



a\E) 



bE + cE 2 , where the constant term a is the 



T taken from [T2], other parameters from [T3] 

Table 4. The parameters r (natural linewidth), Sas (asym- 
metry index), j"° 2 (ratio of intensities between Kai and Ka2) 
and AE (energy difference between Kqi and Ka2) used in the 
fits. 



squared baseline width, i.e. the electronic noise contribu- 
tion. We attribute the b term to statistical fluctuations 
(oc \f~E) in the thermahzation of the AgReC>4 and the c 
term to uncorrected gain instabilities (oc E). 

The fit results for the other parameters are summa- 
rized in Figurel9landfT0l Concerning the Mn peak, a Comp- 
ton back-scattering tail, as discussed in the previous sec- 
tion, accounts almost completely for the observed long 
tail. All other parameters are compatible with the hypoth- 
esis of no energy dependency. The safest conclusion is that, 
if any dependency exists at all, it is hidden by statistical 
and systematic errors or by the cross-correlations between 
the various parameters. 

_ Because of the lower signal-to-background ratio in the 
_RUN15 sum calibration spectrum, a satisfactory fit was 
possible only for the Ti Ka peak. As can be seen in Fig- 
ure[S] and[TUl both the long and short tail parameters are 
found to be compatible with the ones of RUN14. This is 
also illustrated in Figure [TT1 which shows the Ti Ka peaks 
in RUN14 and RUN15, normalized to overlay. 

For what concerns the escape peaks in RUN15, the 
analysis has been restricted to the Re Ka escape peaks 
since the K/3 escape peaks show a too complex structure 
(see FiguresOandU]). The Re Ka escape peaks have been 
"fitted fixing their natural widths r and asymmetry in- 
dexes 5 as, according to the data in Table|U Fits have 
been performed with the two tails as in ([5]) as well as with 
one or no tail, both on the Kal alone and on the two Ka 
peaks together. The results are listed in Tableland shown 
in FigurefTSl for the two Ka peaks. 
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energy [eV] energy [eV] 

Fig. 7. Two examples for fits of the summed spectrum are displayed: Ti Ka (left side) and Mn Ka (right side). The dashed 
lines are the three components in equation ([4]). 



Kai and Ka2 escape peaks 



Peak 
model 


Atot 
[%] 


Ai 
[keV" 1 ] 


^ex P 2 

Atot 

{%] 


A 2 
[keV- 1 ] 


x 2 


no tail 










1.408 


1 tail 






31 ± 14 


35 ±6 


1.224 


2 tails 


1.5 ±0.7 


5.8 ±2.8 


57 ±34 


46 ±6 


1.210 



Table 5. Fit results for the two Kcti and K«2 escape peaks. 



Although the fits with the three models are almost 
undistinguishable by eye, the x 2 suggests that the one 
without tails is the least satisfactory. The two tail and 
one tail model fits have very similar x 2 , although the best 
fit is obtained with the two tail model. 

The results for the Kal escape peak fit are compared 
with the ones for the low energy peaks in Figure^ and 
1101 The comparison confirms that the A parameters do 
not depend strongly on the energy. Considering instead 
the tail amplitudes, it may be concluded that, for the uni- 
formly distributed escape events, whereas the long tail is 
clearly reduced - as expected for a surface effect -, the 
short tail is considerably higher than for the low energy 
peaks. However, it is worth noting that the short tail am- 
plitude determination is hindered by the large escape peak 
width together with the limited statistics, as shown by the 
large error in Table[5] 
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Fig. 8. The fit parameter a is displayed on the left side. On the right side a 2 is displayed in dependence on the X-ray energy. 
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Fig. 9. The fit parameters Ai and A2 are displayed in dependence on the X-ray energy. Filled circles and open squares are the 
RUN 14 and RUN15 data, respectively. 
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are the RUN14 and RUN15 data, respectively. 
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4 Montecarlo Simulations 

To understand the origin of the observed tails another 
Montecarlo simulation based on the Geant4 toolkit was 
developed. Hereby each X-ray line was simulated sepa- 
rately using the energies and line widths listed in Tables 
[2] and |4j respectively. An additional exponential tail was 
added to the left of the Mn lines according to the results 
shown in FigureH] For the simulation of the Re Ka es- 
cape lines, the initial 44 Ti 7-line at E = 78.337 keV was 
used; the escape lines are then automatically generated 
by the Montecarlo simulation. To account for the energy 
resolution of the detector and compare with experimental 
data, the resulting energy spectra were convoluted with 
a Gaussian, whose a was taken from the fitting of the 
respective line. As the AgRe04 crystals do not have a 
regular geometrical shape [51IT7]. they were approximated 
by a cylinder with a diameter of 400/im and a thickness 
of 200^m. 

In addition to the summed spectrum, the sum of the 
two channels 5 and 6 was investigated separately, because 
only these two channels were used for the escape peak 
analysis. 

The Montecarlo simulation describes creation, trans- 
port and interaction of the initial photo-electron created 
by the X-ray photon, of the secondary electrons as well as 
of the secondary X-rays. However, it does not describe the 
energy transfer from the electron(s) to the phonon system 
which ultimately produces the signal observed with ther- 
mal detectors. Incomplete energy detection that creates a 
low energy tail can be caused by either incomplete energy 
transfer from the electrons to the phonon system, or by a 
reduced temperature signal due to phonon energy losses or 
energy deposition in thermally weakly-coupled volumes. 

As the Montecarlo simulation does not account for 
these different processes, the energy Ed ep deposited by 
the electron(s) in each Montecarlo step was multiplied by 
a weight function w(x) < 1.0, where x is the position at 
which the energy deposition happens. 

As the fitting function could be divided into a long and 
a short shoulder, we have used a function w(x) according 
to a model with two loss mechanisms. The first one as- 
sumes that the AgReC>4 crystal has an insensitive surface 
layer of thickness do- The electron energy deposited in 
this layer does not contribute entirely to the thermal sig- 
nal because of an incomplete conversion into phonons. The 
transition from insensitive to sensitive volume is smooth 
and described by a diffusion profile f(d) with a diffusion 
length DL d : 



M = \ 



1 + tanh 



d — dc 



(10) 



The result of this simulation is displayed in Figurc[T3] for 
the case of RUN14. The do and DLd parameters were 



The choice of this function clearly determines the exact tail 
shape. With the present choice the shape will be only approx- 
imately exponential, but the natural linewidth r, the noise a 
and the statistics prevent from appreciating the difference. 



run 


RUN 14 


RUN14 


RUN 15 


spectrum 


sum 


channel 5+6 


channel 5+6 


d [nm] 


20 


52 


52 


DL d [nm] 


64 


52 


52 


r H 


279 


280 


280 


DL r [fj,m] 


1 


1 


1 


L 


0.010 


0.010 


0.010 



Table 6. Optimal parameters for the simulations. The pre- 
cision is Ado = ±5 nm, ADLd = ±2 nm, Aro = ±1 /jm, 
ADL r = ±100 nm, AL = ±0.003. The optimal values and the 
errors were determined by a one-dimensional % 2 analysis. 



found by one-dimensional y 2 analysis to reproduce the 
experimental data in the 1.0 - 7.0 keV energy interval and 
are reported in Table|51 It is shown clearly in Figure fTSI 
that this simple model can describe the long shoulders 
sufficiently well, although it fails to account for the short 
shoulders. 

The Ka peaks in the Montecarlo simulated spectrum 
fFigurerT3"| were fitted with one single exponential tail and 
the resulting long tail amplitudes A exp \ and Ai parameters 
are plotted in Figure[TJ] together with the results shown in 
Figure|9]and[T0l The comparison confirms that this simple 
model indeed reproduces the features of the experimental 
data and, as expected for a surface effect, the long tail 
disappears almost completely in the escape peaks. 

As the short tail parameters A exp 2 and A2 show a com- 
pletely different dependence on the energy, another mech- 
anism must be found to explain the short component of 
the experimental peaks. Two types of effects can be the 
origin of the observed short tail. The first one consists of 
spatial effects, where the tail depends on the energy be- 
cause the spatial distribution of the interaction does. All 
effects which are not related to the interaction position 
belong to the second type. 

In fact there are good reasons to suspect the presence 
of volume effects that are caused by the dependence of the 
thermal signal on the interaction position. The calorime- 
ter is a composite object with an absorber crystal and a 
thermometer attached to it by means of a small glue spot. 
If the phonons do not reach complete thermal equilibrium 
in the whole absorber crystal faster than the detector ther- 
mal time constant, part of them will be lost depending on 
the X-ray impact location^. It is reasonable to assume that 
impact locations farther away from the attachment point 
of the thermometer will create a reduced thermal signal, 
thus creating a local dependence of the energy response. 

Therefore, as an example for this loss mechanism, we 
report here on one model which has a dependence on the 
distance of the energy deposition from the attachment 



4 Indeed the pulses observed with our AgReC>4 mi- 
crocalorimeters cannot be understood only as thermal signals. 
In general they present two decay time constants. The longer 
one - few tens of milliseconds - could be interpreted as the 
detector thermal relaxation to the operating temperature. The 
shorter one - few milliseconds - is probably caused by the di- 
rect interaction of out-of-equilibrium phonons in the attached 
silicon thermometer. 
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Fig. 13. Comparison of the measured (grey) and simulated (black solid line) spectrum for the sum of eight detectors from the 
MIBETA measurement: full spectrum (left side) and Ti Ka calibration peak as a specific example (right side). Whereas the 
long shoulders are simulated sufficiently well, the short shoulders are not accounted for. 




energy [keV] energy [keV] 



Fig. 14. The fit parameters A. £xp \ and Ai are displayed in dependence on the energy of the X-rays. Filled triangles are the 
MC data. 



point of the thermometer. To account for a smooth tran- 
sition of sensitivity, we used a second diffusion model: 



9(r) 



1 



-5 



1 + tanh 



r - r 
DL r 



(11) 



where r stands for the distance of the impact location 
from the attachment point (see FigureQj)] for a schematic 
description), ro is a "thickness" parameter describing the 
gradual decrease of sensitivity for crystal regions farther 
away from the attachment point, and DL r describes the 
transition width. L sets the minimum g value for the r 3> 
ro limit: when L = 1.0, g(r) approaches for r ro. 
The final weighting function w(d,r) is then given by 



w(d,r) = f(d) ■ g(r) 



(12) 



The parameters used in the simulations are summa- 
rized in Table [SJ these parameters correspond to a weight 
function which is 1.0 in the whole crystal volume and de- 
creases just on the lower rim of the cylindrical crystal 
(see FigureQjj]). Figure [16] shows the result of the sim- 
ulated spectrum compared with the measured spectrum 
from Figure [T] 



From the values of the simulation parameters, it can 
also be deduced that the two mechanisms have different 
origin. Whereas do and DLd are rather different for the 
sum of all and the sum of only two channels, this can- 
not be said for the other parameters ro, DL r and L. This 
strengthens our interpretation of the first mechanism as a 
surface effect which is expected to be individual for each 
crystal. The similarity of the other parameters for the dif- 
ferent channels points to a more systematic effect common 
to all crystals. 

Figure [T71 shows the measured as well as the simulated 
Re Ka escape peaks for the 44 Ti measurement (RUN15). 
The Re Ka escape peaks were simulated with the same 
values for ro, DL r and L as for the calibration peaks, 
adding an increased low-energy background from the 44 Ti 
source. The agreement between the experiment and Mon- 
tecarlo as shown in the left panel of Figured?] is fairly 
good. The right panel of FigurefTTl shows the escape peaks 
obtained with the Montecarlo simulation when neither the 
transition natural width nor the detector resolution are in- 
cluded: the plot shows that although no appreciable long 
term tail is expected a very small short term component is 
indeed present in the simulation. However, a comparison 
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Fig. 15. Scheme of coordinates for Montecarlo simulations: d denotes the distance of the event from the detector surface, r 
the distance from the glue point of the thermometer. See text for an explanation of g(r) (equation 1 1 1 p . Units are microns. 



with the results in Table[S] and Figure [TU] shows that the 
Montecarlo fails to predict the correct amplitude of this 
short tail by more than an order of magnitude. 

A model that depends on the distance of the interac- 
tion point from the calorimeter is naturally more depen- 
dent on the geometrical shape of the crystal than a surface 
model. However, other models for spatial dependence with 
different assumptions for the geometry as well as for the 
weighting function g(r) show similar shortcomings. 

The second type of effects which can cause a tail may 
include, for example, complex energy losses due to trap- 
ping in metastable states, event pile-up as well as other 
analysis artefacts. In particular, for what concerns these 
effects, it is quite intriguing to notice that the short tail 
amplitude A exp 2 and the line width <j vary with the energy 
in a similar way. This is illustrated in Figurc[T5] where the 
ratio between the relative amplitude A exv ij A to t and the 
peak broadening in excess of the baseline width {(j\ xc = 
a 2 — <7q) is plotted versus the energy. A common energy 
dependence might be inferred from this plot. This points 
to a possible common explanation for the two peak pa- 
rameters. 

Indeed both the above types of effects could be present 
in our detectors. However, with our current understand- 
ing it is not possible to build models that satisfactory 
explain the short tail in our data. A better understanding 
of the mechanisms causing that short tail would require 
a higher statistics measurement with the 44 Ti source as 
well as calibrations with lower energy (i.e. < 17keV) 7 or 
X-ray sources to better investigate the tail energy - and 
therefore spatial - dependence. 

The important aim of the present investigation was to 
determine to which extent the tails observed in the calibra- 
tion peaks affect the beta spectrum as well. Unfortunately 
nothing conclusive was found for the short tail. Therefore 




10 

energy [keV] 
Fig. 19. Ratio of A exp 2/A to t to a\ 



in order to evaluate the detector response function ex- 
pected for an internal uniformly distributed 2.5 keV beta, 
the Montecarlo was run with few different assumptions 
about the mechanism causing the tail. 

Figure lTSl shows the results. For sake of comparison, 
in the left panel, the Montecarlo simulation was run for 
an external 2.5 keV X-ray with the simple surface effect 
model in (1101) alone (light grey histogram) and with an ad- 
ditional short tail (darker histogram) with the parameters 
found for the CI Kcv X-ray line (see Figurel^l and ITU)) . i.e. 
for instance with an amplitude of about 10%. The simu- 
lations for the internal uniformly distributed 2.5 keV beta 
are shown in the right panel. The simple surface model is 
the light grey histogram at the bottom: the internal beta 
response function still presents a tail due to a small frac- 
tion of beta interactions close to the surface, but its am- 
plitude is about a factor 100 smaller than for the external 
X-ray. The darker histograms correspond to two hypothe- 
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Fig. 16. Comparison of the measured (black) and simulated (grey) spectrum for the sum of eight detectors from the MIBETA 
measurement: the full spectrum on the left side, the Ti calibration peak as a specific example on the right side. For discussion 
see text. 
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Fig. 17. The left panel shows the comparison of the measured (black) and simulated (grey) spectrum for the Re Kai and 
Kq2 spectra. The right panel shows a spectrum without taking into account the Gaussian and natural line widths for a better 
estimation on a possible tail. The light and dark shaded histograms are with and without the inclusion of the surface effect, 
respectively. 




energy [keV] energy [keV] 



Fig. 18. Simulation of the spectrum due to external X-rays (right) and uniform internal ft (left) interactions. The energy is 
2.5 keV in both cases. See text for color explanation. 
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ses for the short tail origin. The black one is caused by 
a volume effect and therefore the tail amplitude is taken 
as found for the Re Ka escape peaks (about 70%). The 
dark grey one is for a non-spatial effect; therefore, it is 
only energy dependent and characterized by the same val- 
ues found as for the external X-rays. In both panels, the 
purely Gaussian peak expected at this energy is plotted 
for reference (solid line). 

The effect of the long tail in the beta response func- 
tion is negligible for an experiment like MIBETA. As far 
as the short tail is concerned, the relevance for MIBETA 
depends of course on which is the correct modeling. This 
lack of knowledge has therefore to be treated as a source 
of systematic uncertainty and correctly quantified. The ef- 
fect of the tails on upcoming higher statistics experiments 
will need a deeper analysis. 



5 Conclusions 
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In this paper we have described a function which allows 
a satisfactory interpolation of the peaks observed in the 
calibration spectra of the MIBETA experiment. 

The results of the fitting show that at least the longer 
of the two exponential tails towards low energy may be 
indeed understood as a result of a surface effect. This 
is confirmed by the Rhenium escape peaks which have 
a strongly reduced long tail. We convincingly tested the 
surface effect hypothesis with a simple Montecarlo simu- 
lation. We used a simple phenomenological model which 
does not pretend to account realistically for the complex 
details of electron and phonon interactions. Nevertheless 
the simulation proves that a surface effect as the one de- 
scribed by (TITJ|) can explain the observed peak shape, no 
matter what is the underlying physics. 

On the other hand, the shorter of the two exponentials 
cannot be understood with a simple assumption. A geo- 
metrical effect of the crystal shape may play a role, but 
for now we can only conclude that we have no convincing 
description of this effect. 

Further investigations with a larger number of regu- 
larly shaped crystals as of in MARE-1 will be mandatory 
to clarify this situation. 
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